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ABSTRACT 


1.  OBJECTIVES 

The  original  objective  for  the  study  was  the  construction  of  analytical  theory  for 
thermal/residual  stresses  in  thin  films  and  bonding  layers.  However,  prompted  by  the 
reviewers  of  the  Journal  of  Applied  Physics,  the  P.I.  involved  one  more  goal,  namely  the 
construction  of  numerical  codes  for  the  2D  and  3D  fields,  which  caused  the  prolongation  of 
the  planned  work  to  the  second  year, 

2.  STATE-OF-THE-ART  IN  THE  FIELD  BEFORE  AND  AFTER  THE  WORK  IS 
DONE 

Thermal/residual  stresses  in  thin  films  and  bonding  layers  are  currently  the  main  problem 
of  reliability  in  high  speed  electronics.  Only  simple  ID  problems  or  heavy  numerical 
experiments  overloaded  by  unnecessary  information  were  available  in  literature  before.  The 
work  was  proposed  to  derive  the  closed  system  of  partial  differential  equations  with  respect 
to  thermal  stresses,  which  would  take  into  account  the  main  simplifying  feature  of  the 
problem,  namely  the  small  thickness  of  films  and  bonding  layers  compared  to  their 
dimensions  in  plan  view.  To  this  end,  the  approach  taken  from  the  dynamics  of  viscous  fluids 
and  pioneered  by  Prandtl  was  used.  As  a  result,  such  systems  and  necessary  boundary 
conditions  were  derived,  providing  the  2D  analytical  theories  for  thermal  elastic  stresses  in 
thin  surface  films  and  thin  bonding  layers.  An  analogous  study  was  performed  for  total 
plastic  deformation  of  thin  films  and  bonding  layers,  and  the  closed  equation  system  of  the 
theory  of  plasticity  was  derived  for  both. 


A  finite  difference  numerical  code  based  on  the  anal5dical  theory  was  developed  and 
verified  in  some  calculations. 

ACCOMPLISHMENTS 

For  the  first  time  the  2D  analytical  theory  of  thermal/residual  stresses  in  thin  films  and 
bonding  layers  was  derived.  This  achievement  is  similar  to  the  construction  of  the  theory  of 
plates  and  shells  by  Kirchhoflf  and  Love  a  century  ago.  The  system  of  partial  differential 
equations  appeared  to  have  a  specific  structure  with  small  parameters  before  the  senior 
derivatives,  which  is  characteristic  for  boundary  layer  problems,  and  preliminary  numerical 
experiments  unveiled  a  boundary  layer  near  the  trim  of  a  thin  film  and/or  bonding  layer  in 
plan  view. 

PERSONNEL 

Luiz  Martinez  took  part  in  the  numerical  calculations  of  thermal  stresses  in  thin  films 
without  a  financial  support  from  this  grant. 

PUBLICATIONS 
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75(2),  pp.  844-849,  1994. 
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any  Peclet  number,  SIAM  J.  on  Applied  Mathematics,  (to  be  published  in  Spring  1997). 
One  section  in  this  paper  treats  an  application  of  the  derived  solution  to  the  growth  of 
planar  voids  in  fine  lines  described  by  the  2D  difl[usion-electromigration  of  vacancies  in 
the  field  of  tensile  stresses  and  temperature. 

6.  NEW  DISCOVERIES  AND  INVENTIONS 

New  effects  of  Poisson’s  ratio  mismatch  and  periodic  law  of  thermal  stress  distribution  in 
plastic  region  were  discovered. 

No  inventions  or  patent  disclosures  were  performed. 
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PARTI 

THE  THEORY  OF  THERMAL  STRESSES  IN  THIN  FILMS 


It  is  generally  recognized  that  thermal  stresses  in  thin  films  are  primarily  responsible  for 
morphological  changes  in  thin  films,  including  hillocks,  whiskers,  and  void  or  pit  formations, 
which  present  serious  problems  of  rehability  in  microelectronics.  In  Part  I,  the  analytical  theory  of 
thermal  stresses  is  constructed  under  some  natural  simplifying  assumptions.  The  latter  include 
that  of  the  thermoelastic  or  thermoelastic-plastic  behavior  of  the  materials  of  a  film  and  substrate 
after  depositing  the  film  and  cooling  the  composite  system  to  an  operational  temperature.  The 
closed  system  of  governing  partial  differential  equations  for  thermal  stresses  in  a  thin  film  is 
derived  for  any  in-plane  shape  of  the  film.  Some  particular  problems  are  solved  in  an  explicit  form 
and  the  implication  of  the  solutions  for  the  prediction  of  hillocks  is  discussed. 


1.  INTRODUCTION 

Suppose,  for  example’,  that  Pb  is  deposited  on  a  Si  substrate  at  room  temperature  in  the 
form  of  a  thin  film  and  then  the  composite  material  is  cooled  to  4.2  K  where  the  Pb  becomes 
superconducting.  The  coefficients  of  the  thermal  expansion  of  Pb  and  Si  are  29.5x10  /°C  and 
2.6xlO'V°C,  respectively.  While  the  Pb  tries  to  shrink,  the  Si  substrate  restricts  it  from  doing  so; 
hence,  in  cooling,  the  Pb  is  under  tension.  These  thermal  stresses  may  include  plastic  yielding  of 
the  film  and  the  formation  of  hillocks.’^"*  An  analogous  phenomenon  is  possible  for  some 
semiconducting  alloys  deposited  on  a  substrate  and  cooled  to  the  nitrogen  temperature  where 
some  semiconductors  become  superconductive.  Upon  heating  the  composite  back  to  room 
temperature,  the  thin  film  tends  to  expand  and  again  is  restricted  by  the  substrate.  Thus,  the  thin 
film  is  under  compression  upon  heating. 


1 


In  Part  I,  the  analytical  theory  of  thermal  stresses  in  thin  films  is  constructed  subject  to  the 
following  assumptions: 

1 .  The  thickness  of  a  film  is  much  smaller  than  the  thickness  of  a  substrate 
and  any  in-plane  dimension  of  the  film  and  substrate. 

2.  No  shear  or  open-mode  cracks  are  formed  during  the  deformation  process 
by  cooling/heating  (ideal  bonding). 

3.  The  material  of  the  substrate  is  homogeneous,  isotropic,  and  linearly 
thermoelastic. 

4.  The  material  of  the  film  is  homogeneous,  isotropic,  and  ideally 
thermoelastic-plastic . 

This  treatment  is  based  on  ideas  advanced  earlier  as  applied  to  simiplier  probieras^ 

2.  SYSTEM  Of  governing  EQUATIONS 
A.  Thin  film 

Consider  a  thin  film  deposited  on  a  substrate  (Fig.  1).  For  clarity,  the  film  is  depicted  in 
Fig.  1  in  plan  view  in  a  rectangular  form  having  dimensions  of  2/xx2/y.  However,  the  theory 
treated  below  is  valid  for  any  film  form.  The  plane,  z=0,  is  chosen  along  the  flat  film-substrate 
interface  and  a  possible  film-unprotected  surface  of  the  substrate.  The  planes,  y=0  and  x=0, 
coincide  with  the  planes  of  symmetry.  It  is  assumed  that  the  boundary  surface  of  the  thin  film  at 
z=t  is  fi'ee  of  tractions  (t  is  the  film  thickness): 

when  z  =  t,\x\<lx,  1^1  <4  cjz  ~  '^zx  -  =  0.  (2.1) 
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when 


The  edges  of  the  film  are  also  fi'ee  of  tractions: 

0<z<tx  =  ±L  =''^xz 

X  X  xy  xz  ^2  2) 

0<z<ty  =  ±ly  Cy  ='C^  =Xyj  =0. 

The  film-substrate  interface  at  z=0  is  ideally  bonded  so  that  all  displacements  and  stresses, 
and  Tzy,  are  continuous  across  the  interface.  Designate  that 


■t^x  ='^2X  Z  =  0;  x^=x^atz  =  0  (2.3) 

In  what  follows,  the  stress,  Cz,  is  considered  negligibly  small  in  both  the  thin  film  and 
substrate  boundary  zone  influenced  by  the  thin  film  and  located  near  the  latter.  Therefore,  the 
equilibrium  equation  with  respect  to  the  z  axis  is  not  essential  and  can  be  omitted. 


Consider  the  equilibrium  equations  with  respect  to  the  x  and  y  axes: 


cx  cy  cz 


CXJ„  cx. 


or, 


^  cx  cz 


=  0. 


(2.4) 


(2.5) 


Here,  Gx,  Oy,  Xxy,  Tzx,  and  %  are  stresses. 

Let  us  integrate  Eqs.  (2.4)  and  (2.5)  over  z  from  z=0  to  z=t.  Using  Eqs.  (2.1)  and  (2.3) 
we  derive 


5a, 


ox 


dy 
da,,  ox 


t  ’ 


X 

y  y 

“  ~  t 


(2.6) 

(2.7) 


oy  ox 


Here,  a^(x,_y),  ay(x,_y),  and  x,^(x,y)  are  corresponding  mean  stresses  in  a  thin  film: 
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FIG.  1.  Schematic  representation  of  a  film-substrate  composite:  (a)  planar  view  of  a  deposited  thin  film,  with 
dimensions,  21^  x  21y  x  t.  (b)  Sectional  view  showing  the  substrate  boundary  zone  with  dimensions  21^  x  21y  x  tb. 


(2.8) 


Prior  to  when  yielding  occurs,  the  film  material  is  thermoelastic  and  satisfies  the  following 


incremental  Hooke’s  law  equations: 


c{At/)  Aa,-u^Ao, 

+  a.A/, 

ax  Ej- 

(2.9) 

c(AF)  Aa^-u^Aa,  _ 

—  +  oc  ^  A  j , 

dy  E,  '  • 

(2.10) 

c(AF) 

dy  ox  Gj- 

(2.11) 
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Here,  U(x,y)  and  V(x,y)  are  mean  displacement  components  along  the  x  and  y  axes,  respectively 
(in  a  thin  film),  E  is  Young’s  modulus,  G  is  the  shear  modulus,  v  is  Poisson’s  ratio,  and  a  is  the 
thermal  expansion  coeflScient.  Subscript  /  corresponds  to  the  film  material.  The  value  of  T 
denotes  the  temperature  with  respect  to  a  certain  reference  temperature  when  all  thermal  stresses 
equal  zero.  (In  the  example  of  the  Introduction,  the  reference  temperature  is  room  temperature, 
at  which  Pb  is  deposited  on  Si.)  The  symbol  A  designates  the  infinitesimal  increment  of  a  quantity 
in  the  thermal  process  of  heating  or  cooling  corresponding  to  the  growth  AT  of  the  temperature. 
The  parameters  a,  E,  G,  and  v,  generally,  are  some  functions  of  temperature.  Equations  (2.6)  and 
(2.7)  hold  for  stress  increments,  too. 

B.  Substrate 

A  thin  film  influences  the  stress  and  strain  field  of  a  substrate  in  a  certain  boundary  zone  in 
the  neighborhood  of  the  thin  film.  Outside  the  boundary  zone,  the  stress  and  strain  field  is  not 
perturbed  by  a  thin  film,  hence  it  is  the  same  field  as  without  the  thin  film,  We  assume  that  the 
boundary  zone  occupies  the  domain  0  >  r  >  -  4  ,  |  x  |  <  4  and  1  |  <  4  >  where  4  is 

the  thickness  of  the  boundary  zone.  The  value  of  4  ,  if  unknown  from  the  physical  formulation  of 
a  problem,  is  considered  as  a  fitting  parameter,  which  can  be  chosen  to  better  approximate  the 
function  of  primary  interest.  For  example,  in  the  case  of  a  rigid,  strip-shaped  thin  film,  it  was 
determined  that^  from  a  comparison  with  the  exact  solution,  4  roughly  equals  one-sixth  of  the 
strip  width. 

Due  to  bonding  with  no  slipping  between  a  thin  film  and  substrate,  and  because  of  the 
small  thickness  of  the  thin  film,  we  can  ignore  the  bending  stresses  in  the  film.  Under  these 
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conditions,  the  thin  film  can  induce  in  substrate  only  shear  stresses,  and  ,  and  in-plane 
displacements,  u  and  v,  and  in-plane  stresses,  ,  Cy  ,  and  %  .  Equilibrium  equations  with 
respect  to  the  x  and  y  axes  have  the  form  of  Eqs.  (2.4)  and  (2.5).  In  the  substrate  boundary  zone, 
we  have 


0  0  d  d 

—  »~  and  — » ~  , 
oz  OK  oz  oy 

so  that  Eqs.  (2.4)  and  (2.5)  take  the  following  form; 

■^  =  0  ,  -^  =  0 
dz  oz 

Also,  in  the  boundary  zone  we  have 


(2.12) 


5u  aw  ov  cw 

~  .  T" » 

oz  OK  OZ  cy 


where  u,  v,  and  w  are  displacement  components,  so  that  the  incremental  Hooke’s  law  can  be 


written  in  the  form 


a(Au)  At, 


a(Av)  At 


yz 


(2.13) 


dz  Gj  ’  dz  Gj 
Subscript  s  corresponds  to  the  substrate  material.  The  A  is  an  infinitesimal  increment  of  a 
quantity  in  the  heating-cooling  process  corresponding  to  the  growth,  AT,  of  the  temperature. 


Integrating  Eqs.  (2.12)  and  (2.13)  over  z,  we  obtain 
='^x  (x,y)  ,  v='^y(^>y)  ’ 

k 

Au  =  z— ^  +  AU(x,y)  ,  Av  =  z— ^  + AV(x,y)  (2.15) 

G,  G, 
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Here,  we  used  the  equilibrium  and  continuity  equations  on  the  film-substrate  interface.  Moreover, 
the  following  compatibility  equations  should  be  satisfied: 

when  z  =  4  Am  =  Amo  ,  Av  =  Avo  .  (2- 16) 

Here,  mo  (x,y)  and  Vo  {x,y)  are  the  corresponding  displacements,  u  and  v,  of  the  field  unperturbed 
by  the  film.  The  functions  mo  and  vo  should  be  determined  fi'om  analysis  of  the  internal  stress  and 
strain  field  in  a  substrate  without  a  film.  In  the  following,  they  are  considered  given.  For 
example,  if  the  substrate  without  a  film  suffers  an  unconstrained  uniform  thermal  deformation,  we 
have 


when  z~  -h  :  Amq  =  -^ct,  AT  ,  Avo  =  AT 
[Here,  we  also  used  the  symmetry  condition  in  Fig.  l.(a).] 
Thus,  from  Eqs.  (2.15)  and  (2.16),  it  follows  that 


AU  =  Au,  + 


AV  =  Av„  + 


G. 


(2.17) 


(2.18) 

(2.19) 


Substituting  AU  and  AV  in  Eqs.  (2.9),  (2.10),  and  (2.11)  by  Eqs.  (2.18)  and  (2.19),  we 


obtain 


^o-x-VfAay  tb  5(AxJ  6{AUo) 

- rr  + 

Ef  Gj  ox  5x 

Aa,-v,.AG,  t,  5(AXy)  a(Avo) 


-ttfAT 


-af  AT 


E.  G3  cy 

f  .  .  ! 

Ax,„  t. 

^  '  7  T  a. 

oy  ox  J  oy  ox 


Gf  Gf 


a(AxJ  5(AXy)l  a(Auo)  g(Avo) 


V 


y 


(2.20) 


(2.21) 

(2.22) 
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a'lAq,  5'  At, 


cy*  oxcy 

G,  ^A^-VfAq 

ttfa  ^  Ef 


g(AVo) 

dy 


-fafAT 


\ 

y 


(2.26) 


Aq  +Aay 

G.  { 


ttb  I  G, 


cx“  cy' 

g(^Uo)^ 

5y  ; 


(2.27) 


The  system  of  three  partial  differential  equations,  Eqs.  (2.25)-(2.27),  will  serve  to  find  the 
increments  of  thermal  stresses  ^ ,  and  f  in  a  thin  film  in  terms  of  x  and  y.  In  view  of  Eqs. 

(2.2)  we  have 
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(2.28) 


x  =  ±lx 

> 

li 

1 

II 

o 

5A  ^  5Ax^y 
dy  ^  dx 

y  =  ±ly 

t> 

II 
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o 

ox  oy 

The  boundary  value  problem  in  Eqs.  (2.25)-(2.28)  is  appropriate  for  both  the  anal5nical 
and  numerical  study  of  thermal  stresses  in  thin  films. 


3.  A  STRIP-SHAPED  THIN  FILM  ON  A  SUBSTRATE 

Consider  the  case  of  Fig.  1(a)  when  4,  »  4  ,  where  a  thin  film  is  deposited  in  the  form  of 
an  infinite  strip,  |  -v  |  <  4  ,  of  width  24  •  In  this  case,  due  to  the  symmetry,  we  have 


—  =0 
5y 


=  0 


(3.1) 


We  also  assume  that  Ai/o  =  ct^xAT  and  Avo  =  ct^yAT.  With  regard  to  Eq.  (3.1),  the 
governing  equation  system,  Eqs.  (2.25)-(2.27),  become 


Here, 


d"  (AcJ  k‘ 

=  TT  A  a,  +  FAT 
dx'  /: 


(3.2) 


k-  = 


ttfiEf 


F  =  -^{l  +  v,.)(af  -aj 


tt 


(3.3) 


(3.4) 


Ac^  =  VcAq,-(af-aJ  E^AT  (3.5) 

The  solution  to  Eq.  (3.2)  meeting  the  boundary  condition,  Aa^  =  0,  when  .v  =  ±  4  ,  is  as 


follows: 
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Aa,  =  FAT-^ 


cosh  (kx  74) 
cosh  k  (T) 


-1 


dT 


(3.6) 


where  F  is  a  constant.  It  is  just  as  easy  to  find  the  solution  for  the  arbitrary  function  F  =  F  (x). 
Since  af ,  a, ,  Ef ,  Vf ,  and  Gs  are  certain  functions  of  temperature  T,  we  integrate  Eq.  (3.6)  over  T 
and  obtain 


5,(T)-q,(T„) 

_  fr  F(T) 
"  %k-(T) 


cosh(xk(T)  /  4) 
cosh  k(T) 


-  1 


dT 


(3.7) 


Here,  To  is  the  initial  temperature  of  the  heating/cooling  process.  Functions,  F(T)  and  k(T)  are 
defined  by  Eqs.  (3.3)  and  (3.4),  e.g.,  in  the  example  problem  of  the  Introduction,  To  is  the  room 
temperature,  T  is  the  helium  temperature,  and 


5>(t,)  =  0  , 

and  from  Eq.  (3.7),  it  follows  that 


•T  F(T) 


a=r- 


■ok-(T) 


cosh(xk  (T)  /  4  ) 
cosh  k  (T) 


-1 


dT 


q=  +  ("la.(T)-a,(T)]E,(T)dT 


+  /. 


=  f 

X  Jj 


•TVf(T)F(T) 
k^(T) 


cosh(xk  (T)  /  4) 
cosh  k  (T) 


-1 


dT 


pF(T)sinh(k(T)x/4) 
■Cx  t ,  k-(T)coshk(T) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 
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Thus,  the  thermal  stresses  are  determined  for  arbitrary  functions,  af  (T),  a,(T),  Ef(T), 


G,(T),  and  v/(T). 


4.  SOLUTION  FOR  LARGE  k 

The  case  of  large  k  defined  by  Eq.  (3.3)  is  of  importance  for  very  thin  films,  which  is 
realized,  for  example,  when  4^  »  ttb  and  Gs  is  of  the  order  of  Ef.  In  this  case,  the  solution  to  the 
governing  equation  system,  Eqs.  (2.25)-(2.27),  constitutes  a  narrow  boundary  layer  that  borders 
the  contour  of  a  thin  film  Everywhere  inside  the  arbitrary  contour  of  a  film,  except  for  the 
boundary  layer,  the  solution  is  given  by  the  following  equations: 


/  = 


-  Vf  Aa  =  EfAT 


0  Au, 

Vox  AT 


-a,. 


AOy  - VfAa^  =  EfAT 


c  Avq 

Vcy  AT 


-a,. 


Ax,„  =  Gf  AT 


^  o  AVf,  o  AUq^ 
^ox  AT  oy  AT 


(4.1) 

(4.2) 


(4.3) 


In  a  normal  case,  when  Auo  =  asXAT  and  Avo  =  cqyAT,  we  have  Aa^  =  Aa^  and  Ax.^^  -  0. 

Inside  the  boundary  layer,  it  is  expedient  to  express  the  thermal  stresses  in  terms  of  Xn  and 
Xt,  where  OxnXt  is  the  local  coordinate  frame  with  the  origin  at  an  arbitrary  point,  O,  of  the 
contour  of  a  film,  and  Xn  and  Xt  are,  respectively,  normal  and  tangential  directions  at  0,  with 
respect  to  the  contour.  In  this  case,  Eqs.  (2.25)-(2.79)  can  be  reduced  to  the  following  equations: 


+  FAT 


(4.4) 
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Ao;,,  =VfAa,^  -(a^  -  ajEfAT 

=  0  ,  Avo  =  a,x,AT  ,  Auo  =  a,x„AT)  ,  (4.5) 


X 


(4.6) 


The  solution  to  these  equations  satisfying  the  boundary  condition,  Aa^^  —  0  when  x  -  0  (on 

the  film  contour)  and  approaching  toward  the  stress  field  given  by  Eqs.  (4.  l)-(4.3),  when  x, 
is  easy  to  find  in  the  following  form: 


[l<(T)/4l' 

=-.f 


-l)dT  , 


'[k(T)/4] 


exp 


x.k(T) 

/ 


dT 


(4.7) 


(4.8) 


The  width  of  the  boundary  layer  is  roughly  equal  to 


k(T)-“^G,(l-v^) 


(4.9) 


While  Xn  is  increasing,  x^n  monotonously  decreases  toward  zero  when  Xn  2c,  and 

monotonously  increases  toward  the  constant  maximum  defined  by  Eqs.  (4.1)  and  (4.2). 
Maximum  values  of  x„  and  follow: 


p  F(T)dT 
k(T)/4- 

fT  F(t)dT 
[k(T)//J' 


(Xn  -^O)) 


Here,  F(T)  and  k(T)  are  defined  by  Eqs.  (3.4)  and  (3.3),  respectively. 


(4.10) 

(4.11) 
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5.  PLASTIC  YIELDING  OF  A  THIN  FILM:  HILLOCK  FORMATION 

A  metal  film  is  elastic  until  it  yields.  Yielding  begins  at  a  certain  temperature  difference,  T 
-  To,  when  the  second  invariant  of  the  stress  deviator  achieves  a  critical  value;® 


q; +C3;,- +3t^ +31:;  =ct;  (5.13) 

Here,  a,  is  the  yielding  stress  under  simple  tension  or  compression. 

According  to  the  solution  of  Section  IV,  the  yielding  can  begin  either  on  the  contour  of  a 

thin  film  or  inside  the  film. (outside  the  boundary  layer).  The  second  case  holds  when  j  in 

Eq.  (4. 10)  achieves  as  first. 


We  consider  plastic  yielding  of  a  thin  film  in  the  most  important  case  of  large  k  at 


d  AUr 


d  Av„ 


■  =  a. 


c  AVq  d  Au^ 


=  0 


(5.2) 


ox  AT  5y  AT  ’  cx  AT  cy  AT 

in  Eqs.  (4.  l)-(4.3),  when  the  following  uniform  biaxial  stress, 

(5.3) 

K  =  o.Vx  =  o)  , 

holds  in  the  thin  film  (except  for  the  boundary  layer)  before  yielding  occurs.  It  is  reasonable  to 
assume  that  the  following  equations: 


q  =  q=a  ,  =  0  ,  =  x  ,  (5.4) 

are  valid  for  the  total  plastic  yielding  of  the  film  (for  sufficiently  large  T  -  To).  Using  the  condition 

of  biaxiality,  Eq.  (5;4),  in  the  yielding  criterion,  Eq.  (5.1),  we  obtain 


cr  +6t-  =Gi 


(5.5) 


According  to  Eqs.  (2.6)  or  (2.7)  the  equation  of  equilibrium  has  the  following  form; 
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do 

t— =  T 

dr 


Here,  r  is  the  radial  distance  from  a  certain  coordinate  origin. 


Let  us  introduce  the  new  function,  g(r) 


o=o,cosg(r)  ,  T  =  -7=03  sin  g(r) 


The  yielding  criterion,  Eq.  (5.5),  is  satisfied  by  Eq.  (5.7),  and  substituting  o  and  x  in  Eq.  (5.6)  by 
Eq.  (5.7)  provides  the  following  equation: 


tV6-^  +  ij  sing(r)  =  0 


This  equation  has  the  following  two  solutions: 

sing(r)=0  g(r)=7m  (n  =0,±1,±2,- •  •) 

0=±0, ,  T  =0 


(5.10) 


Here  C  is  an  arbitrary  constant.  Suppose  that  the  boundary,  r  =  rn  ,  of  a  round  domain,  r  <  th  ,  is 
free  of  tractions: 


when  r  =  rn,  a  =  0 


(5.11) 


Using  Eq.  (5.11)  inEqs.  (5.10),  we  obtain 


TC  ftT 

c=-+-^ 

2  tV6 


and  the  second  solution,  Eq.  (5.10),  becomes 
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(5.12) 


CT=CT,  sm- 


Th  -r 

tV6 


%  -r 


T  — '  COS'" 

V6  tV6 


In  the  process  of  cooling  or  heating,  this  solution,  Eq.  (5.12),  has  developed  from  the 
original  elastic  solution,  Eq.  (5 .3 ),  where  o^  =  ^  =  a  =  ,  x  =  0  at  the  initiation  of  yielding. 


It  is  natural  to  assume  that  the  condition 

cr  =  as,  T  =  0  ,  (5.13) 

is  valid  at  the  center  of  the  round  domain,  r  =  0. 

Using  Eq.  (5.12)  we  find  that  the  condition  of  continuity,  Eq.  (5. 13),  can  be  satisfied  if 

This  is  the  minimal  th  meeting  the  condition,  Eq.  (5.13);  other  rn’s  satisfying  the  same  condition 
are 

r^j  =  TttV3  /  2  +  27mt  Vb  (n  =  1,2,3,- •  •)  (5.15) 

For  large  n,  the  solution  provided  by  Eqs.  (5.12)  and  (5.15)  approaches  asymptotically  to 
the  uniform  solution,  Eq.  (5.9),  in  the  middle  region  that  extends  without  limits,  when  n  is 
increasing. 

The  infinite  series  of  the  solutions  constructed  represent  all  possible  stationary  states  of 
plastic  yielding  equilibrium  of  a  thin  metal  film.  They  probably  provide  the  most  substantiated 
description  of  successive  morphological  changes  occurring  in  thin  metal  films  in  the  process  of 
cooling  or  heating. 


7d;,|-  t  =  3.8476496t 


(5.14) 
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At  first,  the  uniform  solution,  Eq.  (5.9),  is  valid.  However,  the  uniform  yielding  is 
unstable  in  the  post-critical  stage^ . 

The  instability  transforms  the  film  condition  to  the  solution.  Eqs.  (5.12)  and  (5.15),  with  n 
corresponding  to  the  rn  nearest  to  the  original  in-plane  size  of  the  film. 

As  a  matter  of  fact,  the  hillocks  must  initially  have  the  form  of  hexagons,  which  can 
congruently  cover  the  plane  of  a  thin  film  and  maximize  the  released  free  energy*'^.  A  circular 
form  of  hillocks  is  understood  as  a  certain  approximation  to  reality. 

6.  CONCLUSION  TO  PART  1 

The  theory  constructed  opens  a  door  for  predictive  physical  analysis  of  thermal  stresses  in 
thin  films.  The  abundance  of  parameters  significantly  reduces  the  role  of  numerical  experiments 
done  in  this  area  earlier  without  the  analytical  theory. 
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PART  2 

THE  THEORY  OF  THERMAL  STRESSES  IN  THIN  BONDING  LAYERS 


Thin  bonding  layers  made  of  solders  or  metal-filled  adhesives  are  widely  used  in  electronics.  They  are 
mostly  responsible  for  the  integrity  and  reliability  of  computers  and  other  electronic  devices.  Major 
concerns  are  due  to  thermal  stresses  arising  in  the  process  of  electronic  packaging.  In  Part  2,  the 
analytical  theory  of  thermal  stresses  is  constructed  under  some  natural  simplifying  assumptions.  The 
assumptions  include  the  thermoelastic  behavior  of  the  bonded  materials  and  the  thermoelastic-plastic 
behavior  of  a  bonding  layer  after  the  latter  is  deposited  and  the  composite  system  is  cooled  to  an 
operational  temperature.  The  closed  system  of  governing  partial  differential  equations  for  thermal 
stresses  in  a  thin  bonding  layer  is  derived  for  any  in-plane  shape  of  the  layer.  Some  particular  problems 
are  solved  in  an  explicit  form  and  the  implication  of  the  solutions  for  the  prediction  of  planar  voids  and 
cracks  formation  is  discussed. 

1.  INTRODUCTION 

Suppose,  for  example',  that  an  Au-Sn  solder  connecting  an  alumina  ceramic  substrate  and  a  Si 
or  Ga-As  die  (chip),  is  deposited  at  280°C  in  the  form  of  a  thin  layer,  and  then  the  composite  material 
is  cooled  to  room  temperature  of  operation.  The  coefficients  of  the  thermal  expansion  of  the  Au-Sn 
eutectic  alloy,  alumina  and  silicon  are  15.9,  6.7  and  2.8  x  10'^/°C  respectively.  While  the  Au-Sn  tries 
to  shrink,  the  alumina  substrate  and  silicon  die  restrict  it  fi'om  doing  so;  hence,  in  cooling,  the  Au-Sn  is 
under  tension.  These  thermal  stresses  may  induce  plastic  yielding  and  microcracking  of  the  bonding 
layer  producing  a  great  void  or  planar  crack  that  can  split  the  chip  or  otherwise  disintegrate  the 
system^'®. 
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In  Part  2  the  analytical  theory  of  thermal  stresses  in  a  thin  bonding  layer  is  constructed  subject 
to  the  following  assumptions: 

(1)  The  thickness  of  a  layer  is  much  smaller  than  the  thickness  of  a  die  and  any  in-plane 
dimension  of  the  layer  and  bonded  materials. 

(2)  No  shear  or  open-mode  cracks  are  formed  during  the  deformation  process  by 
cooling/heating  (ideal  bonding). 

(3)  The  materials  of  the  substrate  and  die  are  homogeneous,  isotropic,  and  elastic. 

(4)  The  material  of  the  bonding  layer  is  homogenous,  isotropic  and  linearly  thermoelastic- 
ideally  plastic. 

This  treatment  is  based  on  the  approach  advanced  earlier  as  applied  to  simpler  problems^’^ 

The  literature  on  thermal  stresses  is  ample  (see,  for  example.  Refs.  9-19  on  the  subject).  The 
novel  features  of  the  present  approach  compared  with  the  previous  results  of  other  author  as  follows. 

(1)  The  approach  utilizing  the  asymptotic  idea  of  a  boundary  layer  from  hydrodynamics  is  a 
rigorous  and  general  method  of  mathematical  physics.  It  leads  to  the  closed  system  of  governing 
partial  differential  equations  for  all  three  components  of  the  peeling  thermal  stress  in  a  thin  bonding 
layer  in  terms  of  two  in-plane  coordinates  of  the  layer  (Sec.  II). 

(2)  The  equation  system  is  derived  for  any  in-plane  shape  of  the  layer  and  is  formulated  in 
terms  of  stress  increments  as  some  functions  of  temperature.  This  allows  one  to  consider  thermal 
expansions  and  elastic  properties  of  the  constituent  materials  as  some  arbitrary  functions  of 
temperature. 
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(3)  The  solution  to  the  general  equation  system  is  shown  to  generally  split  into  two  solutions- 
one  describing  the  uniform  field  inside  and  in-plane  domain,  and  the  other  describing  the  trim  zone  near 
the  boundary  of  the  domain  (Saint  Venant  effect  zone  or  boundary  layer).  The  first  solution  is  similar 
to  a  potential  inviscid  flow,  and  can  also  be  found  by  elementary  means.  The  second  is  similar  to  the 
viscous  flow  in  a  boundary  layer,  and  can  often  be  found  analytically  for  the  material  parameters  being 
arbitrary  functions  of  temperature  (Secs.  EH  and  EV).  Analytical  results  for  the  Saint  Venant  effect 
zone  in  a  thin  bonding  layer  have  been  unknown  to  this  author. 

(4)  The  well-developed  plastic  yielding  a  thin  bondiog  layer  in  the  internal  domain  is 
considered  using  the  classic  limit  analysis  of  thermal  stresses  (Sec.  V). 

Some  numerical  examples  demonstrating  the  internal  zone  and  boundary  layer  solutions  are 
given  in  Sec.  \T. 

2.  SYSTEM  OF  GOVERNING  EQUATIONS 

A.  Thin  bonding  layer 

Consider  a  thin  bonding  layer  between  a  substrate  and  die  (Fig.  1).  The  plane,  r  =  0,  is  chosen 
along  the  flat  bonding  layer — substrate  interface  and  free  surface  of  the  substrate.  The  x  axis  is  in  a 
cross  section  of  the  bonding  layer  which  can  have  an  arbitrary  form  in  planar  view.  The  x,  y,  and  z 
make  up  a  Cartesian  coordinate  frame. 

The  bonding  layer — substrate  and  bonding  layer — die  interfaces  are  ideally  bonded  so  that  all 
displacements  and  the  stresses,  Cz,  Xr*  and  x^,  are  continuous  across  the  interfaces.  The  following 
designations  are  made; 
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(2.1) 


atz  =  0;  -Zy  =  at  z  =  0 

Xx  =  Xa:  at  z  =  f,  Xy  =  x^  at  z  =  / 

Here  t  is  the  bonding  layer  thickness. 

In  what  follows,  the  variation  of  the  stress,  a^,  is  considered  negligibly  small  in  both  a  thin  bonding 
layer  and  boundary  zones  in  a  die  and  substrate  influenced  by  the  thin  bonding  layer  and  located  near 
the  latter.  Therefore,  the  equilibrium  equation  with  respect  to  the  z  axis  is  not  essential  and  can  be 
omitted. 

Consider  the  equilibrium  equations  with  respect  to  the  x  and  y  axes,  and  let  one  integrate  them  over  z 
fi’om  z  =  0  to  z  =  r.  Using  Eq.  (2. 1)  we  derive 

^CTx  ^^’xy  _  Tx  ~  Xx 

ox  cy  t  ’  (2.2) 

^(Ty  ^  ^Xxy  Xy'Xy 

By  dx  t  ’  (2.3) 

Here,  crx(x,y),cry(x,y),  and  xxy(x,y)  are  corresponding  mean  stresses  in  a  thin  bonding  layer: 

_  1  ^  _  1  ^  _  1  • 

(jx —  ^  Cxdz,  cxy  J  cTydz,  7-,^--  J  rx\'dz 

t  0  ■  ^  0  '  ‘‘o'  (2.4) 


FIG.  1.  Schematic  rq)resentatLon  of  a  substrate — bonding  layer — die  composite  cross  section. 
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Prior  to  when  yielding  occurs,  the  bonding  layer  material  is  thennoeiashc  and  satisfies  the 


following  incremental  Hooke's  law  equations: 

^  ax  E 

(2.5) 

a(AV)^Aay-v'A(^x  +  ^z)^^.j. 
oy  .  E 

(2.6) 

a(AU)  ,  a(AV)  _ 

dy  dx.  G 

(2.7) 

Here,  U  (x,  y)  and  V  (x,  y)  are  mean  displacement  components  along  the  x  and  y  axes, 
respectively  (in  a  thin  bonding  layer),  E  is  Young's  modulus,  G  is  the  shear  modulus,  v  is  Poisson's 
ratio,  and  a  is  the  thermal  expansion  coefficient.  The  subscript  s  will  correspond  to  the  substrate 
material,  subscript  d  to  that  of  the  die,  and  no  subscript  will  correspond  to  the  solder.  The  value  of  T 
denotes  the  temperature  with  respect  to  a  certain  reference  temperature  when  all  thermal  stresses  equal 
zero.  (In  the  example  of  Section  1,  the  reference  temperature  is  the  temperature  of  the  electronic 
packaging  when  the  solder  is  deposited  on  Si,  and  the  chip  is  attached).  Detecting  a  reference  state  is 
often  a  crucial  point  in  the  theory  of  generation  of  residual  stresses^*^.  The  symbol  A  designates  the 
infinitesimal  increment  of  a  quantity  corresponding  to  the  increment  A  T  of  the  temperature.  The 
parameters  a,  E,  G,  and  v,  are,  generally,  some  functions  of  the  temperature.  Equations  (2.2)  and 
(2.3)  also  hold  for  stress  increments. 


B.  Substrate  and  die 


A  thin  bonding  layer  influences  the  stress  and  strain  field  of  a  substrate  and  the  die  in  certain 
boundary  zones  in  the  neighborhood  of  the  bonding  layer.  Outside  the  boundary  zones,  the  stress  and 
strain  field  are  not  perturbed  by  a  thin  bonding  layer;  hence,  they  are  the  same  field  as  with  zero 
thickness  bounding  layer.  We  assume  that  the  boundary  zone  occupies  the  domain,  0  >  z  >  -  4,, 
0  <  X  <  U  and  0  <  >/  <  1^,  inthe  substrate,  and  4  >  r  >  0,  0  <  x  <  1,  and  0  <  >'  <  1^.,  in  the 
die,  in  Fig.  1,  where  4  and  4  are  the  thickness  of  the  boundary  zone  in  the  substrate  and  die, 
respectively,  and  U  3nd  ly  are  some  planar  dimensions  of  the  rectangular  bonding  layer.  The  values  of 
4  and  to,  if  unknown  fi'om  the  physical  formulation  of  a  problem,  are  considered  as  fitting  parameters, 
which  can  be  chosen  to  better  approximate  the  function  of  primary  interest.  For  example,  in  the  case  of 
a  rieid,  strip  shaped  bonding  layer,  it  was  determined^  that  from  a  comparison  with  the  exact  solution,  4 
and  4  are  roughly  equal  one-sixth  of  the  strip  width. 


Due  to  the  small  thickness  of  the  bonding  layer,  we  can  ignore  the  bending  stresses  in  the  layer, 
so  that  the  thin  bonding  layer  can  induce  in  the  substrate  and  die  the  shear  stresses,  Xc  and  Xv:,  and  in¬ 
plane  displacements,  u  and  v,  and  in-plane  stresses,  Ox,  Oy  and  Xxy. 


Using  the  method  of  Section  IIB  of  reference*  we  derive 


AU=.(Au„),+ 
AV  =  (Av„)3+ 


^  =  (Au<,)d  - 

tD^^x 

Gd 

(2.8) 

tD^^'y 

Gd 

(2.9) 
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Here,  (Am<,)„  (Avo)„  (Am^)^,,  and  (Av„);,  are  the  increments  of  corresponding  displacements,  u 
and  V,  in  the  substrate  and  die  of  the  field  unperturbed  by  the  thin  bonding  layer.  These  fimctions  of  x 
and  y  should  be  determined  from  the  analysis  of  the  thermal  stress  and  stress  field  in  the  bonded 
substrate — die  system  with  a  bonding  layer  of  zero  thickness.  In  the  following,  they  are  considered 
given.  The  stress,  cr,  {x,  y),  in  Eqs.  (2.5)  and  (2.6)  should  also  be  determined  from  the  analysis  of  the 
same  system  with  a  bonding  layer  of  zero  thickness. 

From  the  bonding  condition  on  the  substrate — die  interface  with  a  zero-thickness  bonding 
layer,  it  follows  that 


(Aw^)j  (Auo)d  Auo,  (Av^,)^  (Avq)/}  AV(,, 

(2.10) 

and  hence  from  Eqs.  (2.8)  and  (2.9)  we  derive 

4  Gd  a  t.v  =  -  to  Gs  A  T,: 

(2.11) 

4  Gd  AXy  =  -  to  Gs  Axy 

(2.12) 

Writing  Eqs.  (2.2)  and  (2.3)  in  the  incremental  form  and  substitute  At/ 

and  At/  for  Eqs. 

(2. 1 1)  and  (2. 12)  is  obtained; 

1  A  -  _  toGs 

-  A  rx - 

t  toGs  +  tsGD 

1  A  -  _  toGs 

-  A  Tv - ; - 

t  tDGs'^tsGo 


(2.13) 


(2.14) 


Substituting  AU and  AFinEqs.  (2.5),  (2.6)  and  (2.7)  for  Eqs.  (2.8)  and  (2.9),  one  derives 
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(2.15) 


Acrx->'Agy_  vAcr^  ,  t,  3(Auo) 

E  Gs  5x 


+  V  --^-oAT 
5x 


A  o-y  -  vA  o-^  _  y/A  o-z  ^  ts  ^y)  ^  Vo) 

E  E  Gs  5y  5y 


(2.16) 


_  ts 

G  Gs 


c(a  rx)  ^(a  ry) 
5y  5x 


^  g(Auo)  ^  g(Avo) 
dy  ox 


(2.17) 


Substituting  AXv' aii(i  ATv' in  Eqs.  (2.15)  through  (2.17)  by  Eqs.  (2.13  and  (2.14)  provides 


jAaJ  5'(A^y)  toGs  +  t5GE 


8x'  dxdy  tt^  to 


Aax-vAay  vAGj  c(Auo) 


ox 


+  aAT 


,  (2.18) 


(A(Ty)  ^  ^^(Arxy)_  tpGs  +  tsGp 


oy- 


axav 


ttstD 


^Oy-^^Ox  vAcJz  o  (a  Vq)  ^ 


orAT 


cy 


(2.19) 


dxdy 


nv  \ 


-  ^"(Arxy)  ,  _  tpGs  +  tsGp 

ttstD 


/T  + - - -^  + 

^y) 


Ox 


cv 


A  Tx)'  ^Auq)  c(A  Vq) 
G  cy  5x 


(2.20) 


The  system  of  three  partial  differential  equations,  Eqs.  (2.18) — (2.20),  will  serve  to  help  one 

find  the  increments  of  thermal  elastic  stresses  and  %y  in  a  bonding  layer  in  terms  of  x  and  y. 

On  the  boundary  of  a  bonding  domain,  the  corresponding  normal  and  shear  stresses  (tractions)  vanish; 
for  example,  for  a  rectangular  domain  we  have 

x  =  0,  x  =  U  A^  =  Al<y  =  0 

y  =  Q,y  =  \y  A^  =  A1^  =  0  (2.21) 


25 


In  the  limiting  case  when  Go  tends  towards  zero,  the  equation  system,  Eqs.  (2.18)  through 
(2.20),  coincides  with  that  for  a  thin  film*.  In  the  latter  case,  the  ar  is  a  known  normal  load  on  the  thin 
film  surface,  which  is  assumed  to  be  zero  in  ref  8.  The  a.  (x,  y),  Uo  {x,  y,  and  (x,  y)  in  Eqs.  (2. 18) 
through  (2.20)  are  the  values  of  corresponding  stress  and  displacements  on  the  die-substrate  interface 
when  the  thickness  of  the  bonding  layer  equals  zero. 


If  one  considers  the  equilibrium  equation  with  respect  to  the  z  axis  and  integrates  it  over  r  from 
r  =  0  to  r  =  /,  the  application  of  the  Hooke's  law  equation  for  Xc  and  yields 


f 

du  dv 

\ 

—  +  — 

+  t/c 

1 

dx  dy 

J 

i 

^  0 


-  T - 

\cyr  dy~  J 


dz 


(2.22) 


(2.23) 


Here,  u,  v,  and  w  are  the  Cartesian  coordinates  of  the  displacement  vector,  k  is  the  mean  curvature  of 
the  bonding  layer,  and  [A  {x,y,z)]  is  the  discontinuity  of  A  across  the  bonding  layer,  that  is,  [A]  =  A 

{x,y,t)  -  A  {x,y,o).  From  Eqs.  (2.22)  and  (2.23),  it  follows  that  [ctJ  tends  towards  zero  when  t  0, 
because  ii  and  v  are  continuous  along  the  bonded  die — substrate  interface. 


For  small  /,  the  variation  of  a.-  across  the  bonding  layer  is  of  the  order  of  the  negligible 
variation  of  a,  and  ay.  .  Hence,  the  bending  effects  of  the  thin  bonding  layer  can  be  neglected; 
however,  the  bending  of  die  and  substrate,  if  it  is  essential,  is  taken  into  account  by  the  functions,  Uo 
(x,y),  Vo  (x,y)  and  C:  {x,y),  in  Eqs.  (2. 18)  through  (2.20). 
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3.  A  STRIP-SHAPED  THIN  BONDING  LAYER 

Consider  the  case  when  a  bonding  layer  is  deposited  in  the  form  of  an  infinite  stnp,  \x\  <  U,  of 
width  2U  (when  ly  »  U ).  In  this  case,  due  to  the  symmetry,  the  and  are  some  functions  of 
only  X,  and  =  0,  so  that  the  governing  equation  system,  Eqs.  (2.18)  through  (2.20),  become 


Here, 


dH^o-x)^Jc^A-^  +  FAT 

dv  liter  X 


k2  = 


_  literx(l‘ v^)(tDGs'''tsGD) 
ttstoE 


(3.1) 


(3.2) 


k^E 


liter  x(l_  y^) 


a(i  +  v)  ^0^^)  vAAZo 

^  ^  E  AT  ax  AT  -  5y  AT 


(3.3) 


^cry=  K^o-x'^^o'z)'^E.^(Avo)_arEAT  , 


(3.4) 


a(Auo)  ^  g(Avo)_g 
dy  ax 


(3.5) 


To  find  Uo  ,  Vo,  and  a-,  let  us  confine  ourselves  to  a  simple  case  of  a  die  of  rectangular 
parallelepiped  shape  attached  to  a  substrate  of  an  analogous  shape  along  one  face  (Fig.  2).  The 
dimensions  of  the  die  and  substrate  along  the  x  axis  are  assumed  to  be  small  as  compared  to  their 
dimensions  along  the  y  and  z  axes.  An  elementary  analysis  of  thermal  stresses  and  strains  in  this 
composite  system  (with  a  zero  thickness  bonding  layer)  provides  the  following  results: 
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FIG.  2.  A  die-substrate  systems;  (a)  a  cross  section,  and  (b)  a  lateral  view. 


Auo^xAT  aD  v'd(1‘''Fs) 

V^s-V'D^ 

(3.6) 

Avo  =  yATLD+(l  +  v/s)^^=^]  , 

(3.7) 

o 

II 

N 

b 

<1 

(3.8) 

while  the  thermal  stresses  in  the  die  and  substrate  are  as  follows: 

A(J?  =  (1  +  vs)^^^EdAT 

(3.9) 

Acr^  =  (l  +  VD)--— EsAT 

(3.10) 

The  remaining  components  of  the  thermal  stress  tensor  equal  zero. 


The  bending  in  the  Oyz  plane  can  be  neglected  because  the  dimension  of  the  die  and  substrate 
along  the  z  axis  is  of  the  same  or  larger  order  then  the  dimension  along  the  y  axis  (Fig.  2).  In  this 
'f  particular  case,  F in  Eq.  (3 . 1)  is  the  following  constant: 
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The  solution  to  Eq.  (3.1)  meeting  the  boundary  condition,  AQc  —  0,  when  x  —  ±1.*,  after 


integrating  it  over  T,  becomes 


,Tx  =  liter5  /  ^ 

io 


F(T)  cosh(xk(T)xk(T)  /  liter  x) 


k^(T)  I  coshk(T) 


(3.12) 


Here,  To  is  the  reference  temperature  of  electronic  packaging  when  cs:  -  0.  The  F  (T)  and  k 
(T)  are  defined  by  Eqs.  (3.2)  and  (3.11)  where  E,  v,  v„  v^,  a,  a„  a,,,  G,,  and  Go  are  some  known 


functions  of  T.  The  is  defined  by  Eqs.  (3.4,  (3.7),  (3.8),  and  (3.12)  as  follows 


T  a;;  T 

?V  =  I  KT)-^<iT+  I  E, 

To  To 


gs'gp 

Vs  -  V'D  J 


dT 


(3.13) 


The  thermal  stresses  appear  to  substantially  depend,  not  only  on  the  differences  in  the  thermal 
expansion  coefficient  of  die,  substrate  and  bonding  layer  materials,  but  also  on  the  difference  of 
Poisson’s  ratio  of  the  die  and  substrate  materials. 


4.  SOLUTION  FOR  LARGE  k 

The  case  of  large  k  defined  by  Eq.  (3.2)  is  of  importance  for  very  thin  bonding  layers,  which  is 
realized,  for  example,  when  \x  1 4,  and  l/  »  t  4,  and  Gs  and  Go  are  of  the  order  ofE,  In  this  case, 
the  solution  to  the  governing  equation  system,  Eqs.  (2.18)  through  (2.20),  constitutes  a  narrow 
boundary  layer  that  borders  the  contour  of  a  bonding  area.  Everywhere  inside  the  arbitrary  contour  of 
the  area,  except  for  the  boundary  layer,  the  solution  is  given  by  the  following  equations: 
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(4.1) 


A  (Tx  - C7y  =  E  +  —  A  0-2  +  —  ( A  Uo) 


Ao-v_'"Ao-^  =  E 


f  V  0  ^ 

aAT  + — Acr7  +— (Auo) 
~  E  ay  J 


(4.2) 


^5(Auo)  ,  g(Avoy 


A7xv  =  G| 


V  oy 


OX  j 


(4.3) 


In  a  particular  case  of  Section  HI,  described  by  Fig.  2  and  Eqs.  (3.6)  through  (3.10),  we  have 


1 

1 

k: 

1  b 
<3 

(v'-  vd)(1  + Vs)  as-«D 

E 

AT 

l_v 

l-v-  Vs  -  VD 

1 

Acjv 

aD  _  ^  ' 

(l  -  V  vx))(l  +  Vs)  as  -  av 

E 

AT 

l_v 

l-\r  Vs- VD 

(4.4) 


(4.5) 


Inside  the  boundary  layer,  it  is  expedient  to  express  the  thermal  stresses  in  terms  of  x„  and  x,, 
where  Ox„  x,  is  the  local  coordinate  frame  with  the  origin  at  an  arbitrary  point,  0,  of  the  contour  of  a 
bonding  area,  and  x„  and  x,  are,  respectively,  normal  and  tangential  directions  at  0,  with  respect  to  the 
contour.  In  this  case,  Eqs.  (2. 1 8)-(2.20)  can  be  reduced  to  the  following  equations: 


y 

y 


d~(Ao-n) 

dx^ 


at 

liter  n 


A  a,  =  KAJn  +  At7,  )  +  E  _  aEAT 

,  9  liters  (1  -  (tD  Gs  +  U  Gd) 

- -  > 

t  ts  tD  E 


(4.7) 


(4.8) 


(4.9) 
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Fo  = 


k^E 


litern(l-v^) 
5(Aun)  ,  ^Aut) 


a(i  +  vO  +  v)Ao-z  J_Ai^. ^  Aut 
^  ^  E  AT  axt  at  dxn  AT 


^xt  5xn 


-0  ,  Am-- 


tsGo'^tDGs  <^Xn  • 


UsGd 


ATtn“®>  Am- 


-  _  t  to  Gs 


^(^CTn) 


-  rt'  -  0 


(4.10) 

(4.11) 

(4.12) 


’  ‘n  +  dxn 

Here,  subscripts  t  and  n  correspond  to  the  tangential  and  normal  direction  respectively.  The 
Fo^  U  and  k  can  be  considered  to  be  independent  of  .r,  and  Xt.  The  solution  to  these  equations  should 

satisfy  the  boundary  condition,  =  0  when  x„  —  0  (on  the  contour  of  a  bonding  area).  When  x„  — >  oc, 
the  solution  should  approach  the  stress  field  given  by  Eqs.  (4.1)  through  (4.3)  at  the  point  0.  Solving 
Eq.  (4.7)  yields 


o-n 


=  j  Fq  (T)  In  (32  j^g-x„(k(T)//„)  ^  ij  dT 


T.  k"(T) 


(4.13) 


The  Fo  (T)  is  the  value  of  F  at  the  point  0,  The  width  of  the  boundary  layer  is  roughly  equal  to 


^^n(T)  _ UstpE 


k(T)  ^1  -  v^)(tD  Gs  +  ts  Gd) 


(4.14) 


When  x„  is  increasing,  x„*  and  t„‘  monotonously  decrease  toward  zero,  and  c^,  increases  toward  the 
maximum  value 


(-  \  ^  T  FQ(t)  liter  n(T)  . ^ 

'^CTn^max  i 


(4.15) 


Maximum  values  of  interface  shear  stresses,  and  Xn,  follow; 
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T  tsliternFp  Gp 
j^k(tsGD  +  tDGs) 


dT  (xn  =  0) 


LA  =.t  jJDliternFo^dT(xn  =  0) 

V  /max  j  k(tsGD‘^tDGs) 


Here,  Lk  and  Fo  are  defined  by  Eqs.  (4.9)  and  (4. 1 0). 


(4.16) 


(4.17) 


5.  PLASTIC  YIELDING  OF  A  BONDING  LAYER:  FORMATION  OF  PLANAR 
VOIDS  AND  CRACKS 

A  solder  bonding  layer  is  elastic  until  it  yields.  Yielding  begins  at  a  certain  temperature 
difference,  T  -  To^  when  the  second  invariant  of  the  stress  deviator  achieves  a  critical  value"' 


—2  4.— 2  ,  2 - —  —  4-';-2  4.^  2=  2 

cTx  cry  ■^<Jz-crxC7y"crx<7z"cry<^z^'^rxv-  cts 


(5.1) 


where 


2  r"  =  niax 


(^zx)  >  (ck) 


(5.2) 


Here,  Oj  is  the  yielding  stress  under  simple  tension  or  compression.  Yielding  can  begin  either 
on  upper  or  lower  interface  of  a  bonding  layer,  and  either  on  the  contour  of  a  bonding  area  or  inside 


the  latter  (outside  the  boundary  layer).  The  second  case  holds  when  (^)maxr  in  Eq.  (4.15)  achieves 
first. 

We  consider  plastic  yielding  of  a  bonding  layer  in  the  important  particular  case  studied  in 
Section  m  (Fig.  2): 


Designate  that  Qc  =  a. 


Txy  Cz  0,  Xyz  ^ 


(5.3) 
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For  a  well-developed  plastic  yielding,  a  solder  alloy  material  can  be  considered  incompressible; 


we  thus  have 


10 


cry-~(crx'^crz) 


(5.4) 


Additionally,  let  us  assume  that 

-  -  +  -  .r 

'Z'zx  ■'J'zx  ‘  (5.5) 

This  condition  is  particularly  satisfied  if  the  thermal  and  elastic  properties  of  die  and  substrate 

materials  coincide,  and  the  middle  plane  of  a  bonding  layer,  r  =  tH,  is  a  plane  of  symmetiy  of  the 

problem. 


In  this  case,  with  the  account  of  Eqs.  (5.2)  through  (5.5)  the  yielding  criterion,  Eq.  (5.1), 
becomes 

3cr^  +  24r-  =  4(j2  , 

and  the  equation  of  equilibrium  has  the  following  form; 

t— —  =  2r  f-literx  <  X  < +literx) 

dx  ^  (5.7) 

The  ends  of  a  bonding  layer  are  fi'ee  of  tractions  when: 

x  =  ±  U  ,  a=  0  (5.8) 

The  solution  to  the  boundary  value  problem,  Eqs.  (5.6)  through  (5.8)  is  as  follows; 


X  +  liter  X 
tV2 


CTS 

re 


cos 


X  +  liter  X 

xr 


(5.9) 


Here,  21^  should  be  equal  to 


where  «  =  0,  1,  2,  3,....  Particularly,  when  «  =  0  we  have  21;c  iz  t  V2. 

If  21t  <  7r  /  V2,  the  yielding  of  the  entire  bonding  area  is  impossible.  In  this  case,  the  bonding 
layer  is  probably  elastic,  at  least,  partially. 

If  U  satisfies  Eq.  (5.10),  the  stresses,  a  and  x,  in  the  bonding  layer  are  periodical  fimctions  of  x, 
that  is,  they  form  a  firozen  space  wave  of  length,  2  tt  r  V2,  where  the  tension  and  compression  alternate 
during  a  period. 

In  the  middle  of  every  wave  of  tensile  a  the  formation  of  a  planar  rupture,  whose  plane  is 
perpendicular  to  the  x  axis,  is  possible.  The  rupture,  if  it  crosses  interface  can  split  a  die.  It  is 
interesting  that  the  mptures  are  periodically  located  m  the  bonding  layer,  on  the  distance  2  tc  /  V2  fi'om 
one  fi'om  another.  Between  the  nearest  possible  ruptures,  right  in  the  middle,  ct  is  zero  and  t  is 
maximal.  At  these  points  interface  slips  are  possible,  which  are  also  periodic  with  period  k  t  ^2.  The 
ruptures  are  perpendicular  to  the  x  axis  and  can  extend  along  the  y  axis  through  the  entire  bonding 
area. 

If  h  does  not  satisfy  Eq.  (5.10),  in  a  bonding  layer,  there  is  probably  one  segment  of  length  less 
than  K  t  ^2  which  is,  at  least,  partially  elastic.  The  remaining  part  of  the  bonding  layer  is  covered  by  a 
periodic  wave  solution,  Eq.  (5.9). 

The  interpretation  of  an  analogous  periodic  solution  in  Ref  8  with  respect  to  hillock  formation 
in  thin  film.s  is  probably  invalid  in  what  concerns  the  evolution  of  hillocks.  From  the  solution,  it  follows 
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0 

-0.01 

-0.024 

FIG.  3 .  Tlie  sum  of  principal  stresses  a-^+  (surface  graphics). 

that  the  dimension  of  the  hillocks  is  one  and  the  same  from  the  beginning  of  instability  during  the 
process  of  heating  and  cooling.  The  latter  can  only  influence  on  the  level  of  a  plastic  deformation  of  a 
hillock. 

6.  SOME  NUMERICAL  EXAMPLES 

Equations  (2.18)  -  (2.20)  provide  the  closed  system  of  partial  differential  equations  for 
increments  ot  thermal  stresses,  Acr^,  AcTj.  and  A  r.^.,  in  a  thin  bonding  layer  of  an  arbitrarv'  in-plane 

shape.  As  a  numerical  example,  we  consider  the  rectangular  shape  of  the  bonding  layer  in  the  cases  of 
Aar  =  0,  and  Ai/n ,  and  Avd  given  by  Eqs.  (3.6)  and  (3.7),  see  Fig.  2. 
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FIG.  4.  The  (J y.  vs  X  along  y=  0, 


FIG.  5.  The  vs  x  along  y  =\. 

Consider  the  following  input  data: 

gold  -  -  tin  (Au-Sn)  solder 

a=  15.9x  10-^/°C,  v  =  0.35,  G  =  30  G Pa, 

£  =  81  GPa; 

alumina  (AI2O3)  substrate 

a,  =6.7x  10^/°C,  V,  =0.22, 

a  =160  GPa,  £,=390  GPa; 


silicon  (Si)  die 


a£)  =  2.8x  10^/°C,  vd  =0.27, 

GD=150GPa,  £i)  =  381  GPa. 

The  bonding  layer  is  assumed  to  occupy  the  domain,  U  |  <  a  and  |y  I  <  Z),  in  the  xy  plane, 

where 

a  =1.25  mm,  Z)  =  5mm.. 

The  remaining  geometrical  parameters  are  chosen  to  be 
/=100  jam,  ts  =tD  =300  pm. 

The  change  of  temperature,  from  manufacturing  to  operation,  is  taken  to  be 

Ar=200  °c. 

The  author  did  not  have  detailed  data  available  on  the  variation  of  thennal  expansion  and 
elastic  constants  with  temperature  in  this  range,  therefore,  the  material  parameters  were  taken  as 
invariable  in  the  present  numerical  example. 


FIG.  6.  The  ay  vs  y  along  x  =  1. 
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y 


FIG  7.  The  (TyVS  y  along  x  =  0. 

The  calculation  was  done  using  Matematica  and  a  finite  difference  model  with  grid  of  25  x 
100.  The  results  of  the  calculation  are  depicted  in  Figs.  3-7  for  the  quarter,  0  <  x  <  a  and  0  <y  <  b, 
which  is  sufficient  because  of  symmetry. 

The  dimensionless  values 

V  -  _  X  _  y 

E  ^  E  ^  E  a  b 

are  used  in  Figs.  3-7. 

Figure  3  shows  the  three-dimensional  graphics  of  the  sum,  a  ^  +  a  y ,  in  terms  of  x  and  y .  It 

demonstrates  the  uniform  internal  solution  similar  to  a  potential  flow  of  an  inviscid  fluid  and  the 
boundary  layer  solution  around  the  boundary  of  the  rectangular  domain.  Figures  4-7  show  the 
distribution  of  the  characteristic  nonzero  stresses,  (j^and  Gy,  along  the  boundary  contour  and  the 

symmetry  lines  of  the  domain. 
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7. 


CONCXUSION  TO  PART  2 


Thennal  stresses  in  a  bonding  layer  between  a  die  and  substrate  can  be  predicted  and  analyzed 
using  the  equation  system  constructed.  The  analytical  theory  is  necessary  in  this  field  due  to  a  great 
number  of  parameters  in  the  problems  of  such  kind^. 

An  analysis  of  a  particular  problem  discovered  that,  besides  the  difference  of  thermal  expansion 
coefficients,  the  difference  of  Poisson's  ratio  of  the  die  and  substrate  may  also  be  crucial. 

Thermal  stresses  in  a  plastic  state  of  a  bonding  layer  are  shown  to  be  periodic,  which  provides  a 
causative  background  for  possible  formation  of  a  periodic  lattice  of  planar  cracks  or  voids 
perpendicular  to  interfaces. 
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PARTS 

THE  NUMERICAL  CODE  FOR  THERMAL  STRESSES  IN  THIN  FILMS 

AND  BONDING  LAYERS 

Thermal  stresses  are  primarily  responsible  for  moiphological  changes  in  thin  films, 
including  hillocks,  whiskers,  and  void  or  pit  formations,  which  present  serious  problems  of 
reliability  in  microelectronics.  In  Part  3,  a  computerized  code  for  prediction  of  thermal  stresses  in 
thin  film.s  is  suggested.  The  closed  system  of  governing  partial  differential  equations  for  thermal 
stresses  in  a  thin  film  is  utilized  for  any  in-plane  shape  of  the  film.  Some  numerical  experiments 
are  conducted  to  verify  the  code. 

1.  INTRODUCTION 

Suppose,  for  example,  [1],  that  Pb  is  deposited  on  a  Si  substrate  at  room  temperature  in 
the  form  of  a  thin  film,  and  then  the  composite  material  is  cooled  to  4.2K  where  the  Pb  becomes 
superconducting.  The  coefficients  of  the  thermal  expansion  of  Pb  and  Si  are  29.5  x  10  ^°C, 
respectively.  While  the  Pb  tries  to  shrink  the  Si  substrate  restricts  it  from  doing  so;  hence,  in 
cooling,  the  Pb  is  under  tension.  These  thermal  stresses  may  induce  plastic  yielding  of  the  film 
and  the  formation  of  hillocks  [2,  3,  4].  An  analogous  phenomenon  is  possible  for  some 
semiconducting  alloys  deposited  on  a  substrate  and  cooled  to  the  nitrogen  temperature  where 
some  semiconductors  become  superconductive.  Upon  heating  the  composite  back  to  room 
temperature,  the  thin  film  tends  to  expand  and  again  is  restricted  by  the  substrate.  Thus,  the  thin 
film  is  under  compression  upon  heating. 

The  closed  two-dimensional  system  of  governing  partial  differential  equations  for  thermal 
stresses  in  a  thin  film  was  earlier  derived  for  any  in-plane  shape  of  the  film  [5]  based  on  ideas 
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advanced  previously  for  fibers  [6],  The  work  of  other  investigators  in  the  field  on  the  one¬ 
dimensional  theory  of  thermal  stresses  is  acknowledged,  too  [7-13],  Some  general  laws  of 
residual  stresses  and  strains  are  discussed  in  Ref  [14], 

In  Part  3,  the  finite  difference  analogue  model  of  the  governing  equation  system  is 
formulated  for  any  in-plane  shape  of  the  film.  Some  numerical  experiments  are  conducted  for 
rectangular  shapes  which  enable  to  reveal  two  zones,  namely,  the  boundary  layer  zone  on  the  trim 
of  a  film  analogous  to  a  viscous  layer,  and  the  internal  zone  similar  to  a  potential  flow  of  a  fluid. 
The  materials  of  a  substrate  and  film  are  assumed  to  be  homogeneous,  isotropic,  and  linearly 
thermoelastic. 

2.  GOVERNING  EQUATIONS 

Thin  film 

Consider  a  thin  film  deposited  on  a  substrate  (Fig.  1).  For  definiteness,  the  film  is  depicted 
in  Fig.  1  in  planar  view  of  a  rectangular  form  having  dimensions  of  24  x  2/y.  The  plane,  z  =  0,  is 
chosen  along  the  flat  film  substrate  interface  and  a  possible  film-unprotected  surface  of  the 
substrate.  The  planes,  y  =  0  and  x  =  0,  coincide  with  the  planes  of  symmetry.  It  is  assumed  that 
the  boundary  surface  of  the  thin  film  at  z  =  t  is  free  of  tractions  (t  is  the  ;film  thickness): 

when  z  =  t,  |x|  <  4 ,  [yl  <4  ■  •  ^  ^  0  .  (1) 

The  edges  of  the  film  are  also  free  of  tractions: 

when  0  <  -  <  t  x=  ±4  '■  cT:  =  =0 

when  0  <  z  <  t  y  =  ±  ly  :  a :  =  Tyx  =  0.  (2) 
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The  film-substrate  interface  at  r  =  0  is  ideally  bonded  so  that  all  displacements  and  stress  oi, 
and  %  ,  are  continuous  across  the  interface. 

Designate  that; 

Tx=  ts,  at  z  =  0; 

Ty=  at  z  =  0.  (3) 

In  what  follows,  stress  ex,  is  considered  negligibly  small  in  both  the  thin  film  and  substrate 
boundary  zone  influenced  by  the  thin  film  and  located  near  the  latter.  Therefore,  the  equilibrium 
equation  with  respect  to  the  r  axis  is  not  essential  and  can  be  omitted.  Designate  as 
cj^  (x,y),  (7  y  (x,y)  and  (xy)  the  corresponding  mean  stresses  in  a  thin  film: 


y 

t 
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Fig.  1.  Schematic  representation  of  a  film-substrate  composite:  (a)  planar  view  of  a  deposited  thin  film,  with 
dimensions,  24  x  24  x  /;  (b)  sectional  view  showing  the  substrate  boundary  zone  with  dimensions,  24  x  21, b  . 
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Prior  to  when  5delding  occurs,  the  film  material  is  thermoelastic  and  satisfies  the  following 
incremental  Hooke’s  law  equations; 


ox. 


-  Vfts.<yy 


+  a  j-biT 


cibdj) 

— ^ - - - +  a  ftJ 

dy  Ef  ^ 


^AF)  ^AF)  _  ^  xy 
dy  ox.  Gj- 


(5) 


(6) 

(7) 


Here,  U(x,  y)  and  F(x,  y)  are  mean  displacement  components  along  the  x  and  y  axes,  respectively 
(in  a  thin  film),  E  is  Young’s  modulus,  G  is  the  shear  modulus,  v  is  Poisson’s  ratio  and  a  is  the 
thermal  expansion  coefficient.  Subscript,  /,  corresponds  to  the  film  material.  The  value  of  T 
denotes  the  temperature  with  respect  to  a  certain  reference  temperature  when  all  thermal  stresses 
equal  zero.  (In  the  example  of  the  Introduction,  the  reference  temperature  is  room  temperature, 
at  which  Pb  is  deposited  on  Si.)  The  symbol.  A,  designates  the  infinitesimal  increment  of  a 
quantity  in  the  thermal  process  of  heating  or  cooling  corresponding  to  the  growth,  AT,  of  the 
temperature.  The  parameters,  a,  £,  G  and  v,  generally,  are  some  functions  of  temperature.  A 
thin  film  influences  the  stress  and  strain  field  of  a  substrate  in  a  certain  zone  in  the  neighborhood 
of  the  thin  film.  Outside  the  zone,  the  stress  and  field  is  not  perturbed  by  a  thin  film,  hence,  it  is 
the  same  field  as  without  the  thin  film.  This  zone  occupies  the  domain  0  >  z  >  -  tb,  I  ^  1  <  4 
and  !>- 1  <  4 ,  in  Fig.  1(b)  where  h  is  the  thickness  of  the  zone.  The  value  of  tb,  if  unknown  from 
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the  physical  formulation  of  a  problem,  is  considered  as  a  fitting  parameter,  which  can  be  chosen  to 
better  approximate  the  function  of  a  primary  interest.  For  example,  in  the  case  of  a  rigid,  strip 
shaped  thin  film  it  was  determined  that  fi’om  a  comparison  with  the  exact  solution,  4  equals  one- 
sixth  of  the  strip  width  [6],  Because  of  the  small  thickness  of  the  thin  film,  we  can  ignore  the 
bending  stresses  in  the  film.  Under  these  conditions,  the  thin  film  can  induce  in  substrate  only 
shear  stresses,  Xxz  and  Xyz  ,  and  in-plane  displacements,  u  and  v,  and  in-plane  stresses,  ,  CTy  and 

In  the  zone  of  influence  in  the  substrate,  it  was  derived  [5], 

Xxx  =  x,v  (x,  y),  Xy^  =  Xy  (x,  y)  (8) 

Au  =  r^  +  AUf;Y.j;;, 

Av  =  z — —  +  AV(x,y) ,  (9) 

Here,  subscript,  s,  corresponds  to  the  substrate  material;  the  A  is  an  infinitesimal  increment  of  a 
quantity  in  the  heating-cooling  process  corresponding  to  the  growth,  AT,  of  the  temperature. 
Moreover,  the  following  compatibility  equations  should  be  satisfied: 


when  z  =  -tb ; 

Au  =  Auo ,  Av  =  Avo  . 

(10) 

when  z  =  -tb'. 

Auo  =  xas  AT,  Avo  =yas  AT 

(11) 

(Here,  we  also  used  the  symmetry  conditions  in  Fig.  1(a)). 

Thus,  from  Eqs.  (9)  and  (10),  it  follows  that 

AC/  =  A«o+Ml2L  (12) 

Gs 
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tb^r 


y 


AV  —  Avff  + 

The  governing  equation  system  for  T;^has  the  following  form  [5]; 

Aa^-VfA^  d{Auo) 


(13) 


dcdy  ttb 


ox 


dx 


■  +  afAT 


(14) 


0>- 


dx  6y  ttb 


Aa^-vj-Aay  c{Avo) 


0^ 


+  Of  AT 


(15) 


cxdy 


Acr„  +  Act, 


^  c>^{Ar„) 

_  Gs 

Ax^  ^Auq) 

j  a  ’  ^  '’7 

ttb 

- j 

(16) 


The  system  of  three  partial  differential  equations,  Eqs.  (14)  -  (16),  will  serve  to  find  the 
increments  of  thermal  stresses,  0^,0  ^  and  r^y,  in  a  thin  film  in  terms  of  x  and  y.  In  view  of 

Eqs.  (2)  we  have; 


I 


X  =  ±  4 :  Acr.^  =  AXxy  =  0 

y  =  ±ly:  Ady=Arxy  =  0  (17) 

Everyone  of  the  three  equations  (14)-(16)  is  of  the  second  order.  It  means  that  one  more 
^  boundary  condition  equation  is  required  to  close  the  formulation  of  the  boundary  value  problem. 
It  is  reasonable  to  take  the  following  boundary  conditions; 


y 

!1 

1+ 

II 

o 

(18) 

(19) 
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Physically  it  means  that  the  shear  force  on  the  film  border  tangential  to  the  film  in-plane  contour  is 
neglected.  The  boundary  value  problem  in  Eqs.  (14)-(19)  will  be  used  for  the  numerical  study  of 
thermal  stresses  in  thin  films. 

3.  FINITE  DIFFERENCE  ANALOGUE  OF  THE  BOUNDARY  VALUE  PROBLEM 
Introduce  the  following  notation: 

^  =  ;7  =  y  (20) 


X= - a^aflSl 


d{^UQ) 


h^ajbl- 


^  ^Avq) 
dy  8x 


The  governing  equation  system  (14)  -  (16)  in  this  notation  can  be  written  as: 


d^h  . 
+  -  =/l: 

d^dn 

d^arj 

a^h  a'^'ci 

-  2 
CT] 

(g-Vff  +  b) 
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The  boundary  condition  equations  for  the  rectangular  domain  I  ^  I  <  ^  I  plane 

follow  from  Eqs.  (17)-(19): 

#  =  ^4:  =  |-|  =  0  (26) 

n  =  ±ly.  g~h^o.  f+*=0  x(4  =  /,//,/,  =  /v/')-  (22) 

Let  us  choose  the  finite  difference  rectangular  grid  with  grid  cell  size,  and  At],  measured  in 
whole  numbers.  Designate: 

fn,.n  =f(m,n),  g„,„  =g(m,n),  h„_„  =  h(m,n).  (28) 


FIG.  2  The  3-D  graphics  of  the  sum  .  in  terms  of  x  and  y . 
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I 


FIG.  3.  The  maximum  shear  stress  I  ct,  -  ct,  +  In^y  I  /2,  in  terms  of  x  and  y . 

The  value  of  and  Ari  less  than  1  has  no  meaning,  because  the  governing  equation  system  (14)- 
(16)  loses  sense  when  the  characteristic  size  of  a  domain  in  the  xy  plane  has  an  order  of  the  film 
thickness. 

In  the  finite  difference  form  the  equation  system  (14)-(16),  by  taking  Aq  =  1  and  Ar|  =1, 
can  be  written  as  follows: 


fm+ln  -  '^fm.n  f m-\.n  +  ^m.n  +~^/w+l./7-l  ^  ^tn^\.n  -  ^m-ln  ^  ^m.n+l 


hm.n-\  =  ^[fm.r,  “  V^m.n  +  «)  i 


(29) 


Sm.n+\  2^7n.  W  +  ^  ^  ^  ^m+\.n  2 


“  ~  ^m.n+\  ~  ^mn-\  ~  '^f  fm.n  +  • 

1  1  1  1  * 

-  — (/  +  ^)ot,/j+1  ^  ~  ^^m.n  +  ^m-\.n  +  ^m.n+l  ~ 


(30) 


+  ^m.n-\  -  ^^[2(1  +  Vy  )h^  „  -  c  . 
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(31) 

The  boundary  condition  equations,  (26)  and  (27),  are; 

m^M:  fun  =hun  ^  0,  -  gU.n  =  hu-Vn  -  hu.n  (32) 

n=N:  gm!^  =/»m/V  =  0,  fm+UV-fmi/  =  -hmN  ^  h„^ .  \  .  (33) 

4.  NUMERICAL  EXAMPLE 

To  examine  the  numerical  model  (29)-(33)  let  us  consider  the  normal  isothermal  expansion  with 
no  additional  loadings.  In  this  case: 

a  =  b( a  j- -  a ^)ls.T,  c  =  0  .  (34) 

Consider  the  following  data; 
silicon  substrate  (Si) 

a.  =  2.6  X  IQ-^'C;  v,  =  0.27;  G.  =105  Gpa; 
lead  thin  film  (Pb) 

a/  =  29.5  X  lO-^'C;  v,  =  0.37;  Ef  =  80  Gpa; 
film  geometry  (25  x  100  rectangular  domain) 

ly=N  =  5Q-,  t  =  \  /mi;  =  3  /rni; 

The  temperature  difference,  from  manufacturing  to  operation,  is 

AT=  300°C 

Hence,  we  have 

A=-^  =  0.625,  fl=^  =  8xl0■^ 
c  =  0  Vy  =  0.37, 

The  authors  did  not  have  available  detailed  data  on  the  variation  of  thermal  expansion  and  elastic 
constants  with  temperature  in  this  range,  therefore,  the  material  parameters  were  taken  invariable 
in  the  present  numerical  example. 


51 


aoom 


0.0006 


0.0002 


FIG.  4.  Stress  vs.  x  for  =  1,  0  <  x  <  1. 


FIG.  5.  Stress  vs.  x  for  >'  =  0,  0  <  x  <  1. 


The  calculation  was  done  using  Matematica  and  a  finite  difiference  model  with  grid  25  x 
100.  The  results  of  the  calculation  are  depicted  on  Figs.  2-7  for  the  quarter,  0  <  x  <  a  and  0  <>'  < 
b,  which  is  sufficient  because  of  symmetry.  The  dimensionless  values 
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(35) 
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-0.0124 

FIG.  6.  Stress  a*  ^  vs  y  fox  x  =0.0  <y  <\. 

-0.0036 

-0.0015 


FIG.  7.  Stress  vs.  y  for  x  =  1,  0  <  >'  <1. 

were  used  in  Figs.  2-7. 

Figures  2  and  3  show  the  three-dimensional  graphics  of  the  s  um,  Ox  +  cjy ,  and  maximum 
shear  stress  in  terms  of  x  and  y  .  They  demonstrate  the  uniform  internal  solution  similar  to  a 
potential  flow  of  an  inviscid  fluid  and  the  boundary  layer  solution  around  the  boundary  of  the 
rectangular  domain.  Figures  4-7  show  the  distribution  of  the  characteristic  non-zero  stresses,  Gx 
and  Gy  ,  along  the  boundary  contour  and  the  symmetry  lines  of  the  domain. 
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5.  CONCLUSION  TO  PART  3 

A  numerical  code  was  suggested  to  find  thermal  stresses  in  thin  films  and  bonding  layers.  An 

examination  of  the  model  proved  that  a  narrow  boundary  layer  exists  in  the  neighborhood  of  the 

border  of  the  domain  occupied  by  the  thin  film  or  bonding  layers. 
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